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We consider the non-equilibrium dynamics of disordered systems as defined by a master equation 
involving transition rates between configurations (detailed balance is not assumed). To compute the 
important dynamical time scales in finite-size systems without simulating the actual time evolution 

■ which can be extremely slow, we propose to focus on first-passage times that satisfy 'backward master 
equations'. Upon the iterative elimination of configurations, we obtain the exact renormalization 
rules that can be followed numerically. To test this approach, we study the statistics of some first- 
passage times for two disordered models : (i) for the random walk in a two-dimensional self-affine 
random potential of Hurst exponent H, we focus on the first exit time from a square of size L X L if 
one starts at the square center, (ii) for the dynamics of the ferromagnetic Sherrington-Kirkpatrick 
model of N spins, we consider the first passage time tf to zero- magnetization when starting from a 

^— i ' fully magnetized configuration. Besides the expected linear growth of the averaged barrier kit/ ~ N, 

, we find that the rescaled distribution of the barrier (ln£/) decays as e - "'' for large u with a tail 

■ exponent of order n ~ 1.72. This value can be simply interpreted in terms of rare events if the 
sample-to-sample fluctuation exponent for the barrier is 4>width = 1/3. 
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INTRODUCTION 



In statistical physics, any large-scale universal behavior is expected to come from some underlying renormalization 
('RG') procedure that eliminates all the details of microscopic models. For the non-equilibrium dynamics of disordered 
systems, we have recently proposed a strong disorder renormalization procedure in configuration space that can be 
defined for any master equation [H, 0, Q : it is based on the iterative elimination of the smallest barrier remaining 
in the system, and thus generalizes the real-space strong disorder procedures that had been previously defined for 
random walks in one-dimensional random media [2, 0, 0, 0, HI • However, as for all strong disorder renormalization 
procedures (see |9j for a review), the results are asymptotically exact only near "Infinite disorder fixed points" 
: for the dynamical problems defined by a master equation, this means that the strong disorder renormalization 
procedure will give asymptotically exact results only if the renormalized distribution of barriers becomes broader 
I , and broader upon iteration (see [l| for a more detailed discussion). In the present paper, we show that one can 
^ ■ obtain exact renormalization rules, without any strong disorder hypothesis, if one considers the 'backwards master 
0^ \ equation' satisfied by first-passage times. It turns out that the renormalization rules for the transition rates are 
formally identical to the strong disorder rules introduced in [H, 0, but the interpretation, the goals, and the validity 
of the two approaches are different, as we explain in more details below. 

From a numerical point of view, the main limitation of Monte-Carlo dynamical simulations of disordered systems is 
that the dynamics in the presence of quenched disorder becomes extremely slow as the system size increases (see for 
instance the introduction of our recent work [Io| and references therein) . It is thus important to develop other methods 
to characterize the dynamical properties of disordered systems without simulating the dynamics. For instance in our 
previous work [Toj , we have proposed to use the mapping between any master equation satisfying detailed balance 
and some Schrodinger equation in configuration space, to obtain the largest relaxation time of the dynamics via 
any eigenvalue method able to compute the energy of the first excited state of the associated quantum Hamiltonian. 
Here we propose another strategy based on the 'backwards master equation' satisfied by first-passage times. The 
fact that first-passage times satisfy 'backwards master equation' is of course very well-known and can be found in 
most textbooks on stochastic processes (see for instance [III EH [H, E3])- In the field of disordered systems, the 
backward Fokker-Planck equation has been very much used to characterize the dynamics of a single particle in a 
random medium (see for instance [!, EE EE E3, EE EH ) , but to the best of our knowledge, this approach has not yet 
been used in higher dimension, nor for many-body problems. To test the present approach, we compute the statistics 
of first-passage times over the disordered samples of a given size for two disordered models (i) a random walk in a 
two-dimensional random potential (ii) a mean-field spin model. 

The paper is organized as follows. In section[Hl we recall that first-passage times satisfy 'backward master equation'. 
In section IIII1 we derive the corresponding renormalization rules and discuss the similarities and differences with 
respect to strong disorder renormalization procedures. We then apply this approach to two types of disordered 
models : section IIVI concerns the problem of a random walk in a two-dimensional self-affine potential, and section [V] 
is devoted to the dynamics of the ferromagnetic Sherrington-Kirkpatrick model. Our conclusions are summarized in 
section IVH 
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II. REMINDER ON FIRST-PASSAGE TIMES AND BACKWARD MASTER EQUATIONS 



A. Master equation denning the stochastic dynamics 



In statistical physics, it is convenient to consider continuous-time stochastic dynamics defined by a 'forward' master 
equation of the form 

= (C) W (C - C) - P, (C) W out (C) (1) 

that describes the evolution of the probability Pt(C) to be in configuration C at time t. The notation W (C — > C) 
represents the transition rate per unit time from configuration C to C, and 

W out (C) = Y,W(C^ C) (2) 

a 

represents the total exit rate out of configuration C. 



B. Backward Master Equation satisfied by the first-passage time 



Let us now focus on the following problem : suppose the dynamics starts at t = in configuration C, and one is 
interested in the random time t where the dynamics will reach for the first time any configuration belonging to a 
given set A of 'target' configurations. As is well known (see for instance the textbooks [Til [l2l [T3L IT3|). the mean 
first-passage time t^ a \C) =< t > (where the notation < . > represents the average with respect to the dynamical 
trajectories) satisfies the following 'backward master equation' for all configurations C not in the set A 

^ W (C - C ') t{A) ( C ') - W »ut (C) (C) = -1 (3) 
c ' 

whereas all configurations in the set A satisfy the boundary conditions 

tW(C6A)=0 (4) 

The derivation of Eq. [3] consists in considering what happens during the first time interval [0, dt] if the system is in 
configuration C at t — : at time dt, the system is either in configuration C with probability [W (C —* C) dt], in which 
case the remaining mean time is t^ a \C '), or the system is still in configuration C with probability [1 — W ou t (C) dt], 
in which case the remaining mean-time is (C ) . By consistency, the mean first passage time has thus to satisfy at 
first order in dt 

(C) =dt + J2[W(C - C) dt] (C ') + [!- W ou t (C) dt] (C) (5) 
c ' 

yielding Eq. [3] 

The backward master equations of Eq. [3] can be solved numerically by any method appropriate for linear equations 
with fixed right hand-side. In the next section, we show that they satisfy exact renormalization rules. 



III. RENORMALIZATION RULES FOR FIRST-PASSAGE TIME PROPERTIES 



A. Iterative elimination of configurations 



If one eliminates iteratively the configurations from the system of Eqs [3J satisfied by the first-passage times, the 
renormalized equations for the surviving configurations keep the same form, but with renormalized transition rates 
W and renormalized right-hand sides K R 

J2W R (C - C) T M (C ') - W Q R ut (C) T M (C) = -K R (C) (6) 
c ' 



3 



This equation for C — Co can be used to eliminate (Co) via 

1 



(Co) 



W R ut (C ) 



^W R {Co -^C")t {A) {C ")+K r {Co) 



(7) 



Upon the elimination of the configuration Co, the renormalized coefficients W R and K R evolve according to the 
following rcnormalization rules for the surviving configurations C 



w 



Rnew 



W£t (Co) 
(Co) 

K «^ (c) =jrfl(c) + w^ ^°V (c ) 



(8) 



B. Renormalization rules for other observables satisfying 'backward master equation' 

Since other observables are known to satisfy similar 'backward master equations', it is interesting to discuss here 
their renormalization rules and to compare with Eqs [5J 

1. Higher moments of first-passage times 

Above we have considered the first moment t^(C) =< t > of the first-passage time in the set A when starting in 
configuration C. However, one may consider the higher moments (C) =< t n > that satisfy the following 'backward 
master equation' [ll|, [l2|, [l4[ ) for all configurations C not in the set A 

J2 W (C - C) r n (C ') - W out (C) r n (C) = -nr n _j (C) (9) 

C ' 

whereas all configurations in the set A satisfy the boundary conditions 

r„(C e A) = o (10) 

The derivation of Eq. [5]consists again in considering what happens during the first time interval [0, dt] (see explanations 
before Eq[5]). The higher moments of first-passage times can be thus computed one after the other : if one knows the 
moments of order (n — 1), one can compute the moments of order n via the same renormalization rules of Eq. [5] : 
the only change will be in the initial condition for the right handside that will read K n niUal (C) = nr n ^i (C) instead 
of K™f al (C) = 1. 



2. Escape probabilities 

The simplest quantities that satisfy some backward master equation are the escape probabilities. Suppose the 
dynamics starts at t = in configuration C, and one is interested into the probability Eb/a{C) to reach first any 
configuration belonging to a set B of configurations before any configuration belonging to another set A of configura- 
tions. As is well known (see for instance the textbooks [lli E2, 0, 0] ) , this escape probability E B / A (C) satisfies the 
following 'backward master equation' for all configurations C neither in the set A nor in the set B 

W (C -> C) E B/A (C ') - W out (C) E B/A (C) = (11) 

c ' 

whereas the configurations in the set A or in the set B satisfy the boundary conditions 

E B/A (C g A) = (12) 
E B/A (C G B) = 1 (13) 

(14) 



The backward master Eq. [TT] does not contain any right handside in contrast to Eq. [3] : the iterative elimination 
of configurations will lead to renormalized transition rates that follows the same two first rules of Eq. [5] 
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C. Similarities and differences with the strong disorder renormalization of Refs 

It turns out that the renormalization rules for the transition rates given in the two first lines of Eq. [3] are formally 
identical to the strong disorder rules introduced in [3, 0] . It is thus important to stress here why the interpretation, 
the goals, and the validity of the two approaches are significantly different : 

(i) The present renormalization rules are exact for any dynamics defined by a master equation. But they yield results 
only for observables like first-passage times that satisfy backwards master equations with fixed right hand-side. 

(ii) On the contrary, the strong disorder renormalization procedure introduced in [l], [3] aims to renormalize the 
forward master equation of Eq. [TJ i.e. the full time evolution of the probability distribution Pt{C). It will become 
asymptotically exact at large times only for dynamics governed by an 'infinite disorder fixed point' (see more details 
in However whenever it is the case, it can yield results for any universal observable (i.e. exponents or rescaled 
distributions). 

IV. RANDOM WALK IN A TWO-DIMENSIONAL SELF-AFFINE POTENTIAL 
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FIG. 1: (Color on line) Statistics of the first exit time t ex u from a square of size L x L when starting at the center for the 
random walk in a self-amne random potential of Hurst exponent H — 0.5 : (a) Probability distribution QL^Texit ~ In t ex u) for 
L — 20, 30, 40, 50, 60, 70 ; (b) the log- log plot of the disorder-average T e xit(L) = lnt ex it(L) corresponds to the barrier exponent 
ip = H = 0.5 (Eq. HOJ). 

In this section, we apply the method of the previous section to the continuous-time random walk of a particle in a 
two-dimensional sclf-affine quenched random potential of Hurst exponent H = 0.5. Since we have studied recently in 
Q the very same model via some strong disorder renormalization procedure, we refer the reader to Q and references 
therein for a detailed presentation of the model and of the numerical method to generate the random potential. Here 
we simply recall what is necessary for the present approach. 

We consider a two-dimensional square lattice of size L x L. The continuous-time random walk in the random 
potential U(r) is defined by the master equation 



dP t (f) 

dt 



= J2 P t ( f ') W{f'^r)- P t (?) W out (r) (15) 



where the transition rates are given by the Metropolis choice at temperature T (the numerical data presented below 
correspond to T = 1) 

W(r^r') = 5 <r ~ P r > min(l,e-^')^»/ T ) (16) 



where the factor <5 < f i ^ /> means that the two positions are neighbors on the two-dimensional lattice. The random 
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potential U(r) is self-affine with Hurst exponent H = 0.5 



[U{r)-U{f')f 



<i\2H 



(17) 



We focus here on the first-passage time t^ a \Cq) corresponding to the following conditions : (i) the initial con- 
figuration Co is the center of the square (xq — L/2,yo — L/2) (ii) the set A of 'target configurations' is the set of 
all boundary sites of the square, i.e. having x = 1, x = L, y = 1 or y = L. The first-passage time t^ a \Cq) thus 
corresponds here to the first exit time t ex u from the square L x L when starting at the center. The appropriate 
variable is actually the barrier defined as 



In t„ 



(18) 



On Fig. [T](a), we show the corresponding probability distribution Ql^cxU = ln^eait) for various sizes 20 < L < 80 
with a statistics of 9.10 5 > n s (L) > 36. 10 2 disordered samples. 

As shown by the log- log plot of Fig. [T](b), we find that the disorder-averaged value T exit (L) scales as 



with a barrier exponent tp of order 



T exlt {L) cx 

L — >oo 



ip = H = 0.5 



(19) 



(20) 



These results are in agreement with scaling arguments on barriers |2Gi . [2l| , with the strong disorder renormalization 
approach of [|[ , and with the computation of the relaxation time to equilibrium [l(| • 

V. DYNAMICS OF FERROMAGNETIC SHERRINGTON-KIRKPATRICK MODEL 
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FIG. 2: (Color on line) Statistics of the first time tfu p where the magnetization vanishes, for the ferromagnetic Sherrington- 
Kirkpatrick model of N spins (Eg [22)) : (a) the disorder- average la tfu p (N) grows linearly with N (Eq. |25|| . The disorder-average 
lnt eq (N) associated to the largest relaxation time t eq (N) towards equilibrium as computed from the method of Ref [Hi] is also 
shown for comparison, (b) The rescaled probability distribution Qfu v [u) of Eq. I27I shown here in log scale to see the tail of 
Eq. I28I exactly coincides with the rescaled probability distribution Q eq (u) as computed from the method of Ref [TTJ] : the tail 
exponent is for both of order rj ~ 1.72 (Eq. 1291 ). 

As an example of application to a many-body disordered system, we consider in this section the ferromagnetic 
Shcrrington-Kirkpatrick model where a configuration C = {Si} of N spins Si = ±1 has for energy 

U = - J v S i S i ( 21 ) 

l<i<j<N 
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where the coupling Jy between two spins Si and Sj contains a non-random ferromagnetic part Jo and a random 
Gaussian part J^ of zero-mean Jy = and variance unity Jy = 1 with the appropriate mean-field rescalings 

si he mm 



j = + J ^ (22) 



Here we consider the values Jq = 2 and temperature T = 1 where the model is in its ferromagnetic phase [22l . |23| . 
H3, [25| to study its dynamical properties. The Metropolis dynamics corresponds to the master equation of Eq. Q] in 
configuration space with the transition rates 

W (C - C) - 6 <c ,c ■> min (l, e"W c ')- u WV T ') (23) 

where the factor <5<c,e '> means that the two configurations are related by a single spin flip. 

We focus here on the first-passage time t^ a \Cq) corresponding to the following conditions : (i) the initial configu- 
ration Co is the fully ferromagnetic configuration of magnetization Mn = Xh=i Si = N where all spins arc Si = +1. 
(ii) the set A of 'target configurations' is the set of all configurations of zero magnetization Mn = YliLi Si = 
(we consider only even N). The first-passage time t^ a \Cq) thus corresponds here to the first time tfu p where the 
magnetization Mjy vanishes. 

We have computed the distribution Ql of the barrier defined as 

r 'flip = lni/Kp (24) 

over the disordered samples of even sizes 4 < N < 12 with a statistics of 2.10 s > n s {L) > 6.10 2 samples. As a 
comparison, we have also computed the distribution of the barrier T eq = \nt eq , where t eq is defined as the largest 
relaxation time towards equilibrium via the method described in our previous work [10| . Since the system is in its 
ferromagnetic phase, one expects that the disorder-average of the barrier grows as 



Tfnp(N) = hit fup cx N (25) 

TV— >oc 



and this is indeed what we measure both for Tfu p (N) and for T eq (N) as shown on Fig. [2] (a). The width A(iV) of 
the barrier distribution is expected to grow with a subleading exponent < ip W idth < 1 

AW = (r% P (N) - (r>-(7V)) 2 ) 1/2 ^ N*-«» (26) 

but we are not aware of any theoretical prediction or any previous numerical measure of this sample-to-sample 
fluctuation exponent ip 'width- This is in contrast with the spin-glass Sherrington-Kirkpatrick model correspondin g to 
Jo = 0, where the barrier exponent has been much studied either theoretically [2(1 [2?| or numerically [Tol. T28L l29l l3dl. 

MM- 

With our numerical data limited to small sizes 4 < N < 12, we see already the expected linear behavior of the 
disorder- average of Eq. [2H]as shown on Fig. [5] (a), but we are unfortunately not able to measure the exponent i/j W idth 
of Eq. [55] from the variance. However, since for these small sizes we can study a large statistics of disordered samples, 
we have measured the rescaled distribution Q defined as 



QdTfUp) - -^QfUv [u = ^ J (27) 

We find that the rescaled distribution Q(u) shown on Fig. [5] (b) presents at large argument the exponential decay 

lnQ fUp (u) cc -vP (28) 

with a tail exponent of order 

r] ~ 1.72 (29) 

We have moreover checked that the rescaled distribution Qfu p {u) exactly coincides with the rescaled probability 
distribution Q eq (u) as computed from the method of Ref [lfj. 
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To interpret the value of Eq. [29] one may propose the following rare-event argument. Since the system is in its 
ferromagnetic phase, it seems natural to expect that the anomalously large barriers in the dynamics will correspond to 
the samples that have anomalously strong ferromagnetic contributions coming from the random parts of the couplings 
in Eq. [22]: with an exponentially rare probability of order e~( cst ' N , the N 2 random variables Jij will be all positive. 
Then instead of being finite, the local field hi = J2j JijSj on spin Si will be of order TV 1 / 2 , and one thus expects a 

barrier of order N 3 / 2 . If one plugs these values in Eqs 1271 and l28l one obtains, for the powers of N in the exponentials, 
the consistency equation 



3 
2 

For instance ipwidth = 1/2 would correspond to r\ — 2. The value 



ipwidth )v = 2 (30) 



ipwidth = \ (31) 



would correspond to the tail exponent value 



vtywWh = ~) = y = 1-714... (32) 

which is extremely close to the value that we measure numerically (Eq. |2"9"]) . A tentative conclusion would thus be 
the following : at the small sizes that we can study, we cannot measure the width exponent ipwidth from the variance, 
but we can measure the tail exponent rj that contains the information on ipwidth if one can properly identify the rare 
events that dominate the tail. In the ferromagnetic phase considered here, we believe that the rare events dominating 
the tail are the anomalously strong ferromagnetic samples described above, so that our measure of the tail exponent 
of Eq. [29] would point towards the value of Eq. [31] for the width exponent. Of course, this type of indirect reasoning 
based on rare events remains rather speculative, and a direct measure of ip W i d th from the variance for large sizes N via 
Monte-Carlo simulations would be very welcome (to the best of our knowledge, the variance has only been measured 
up to now for the case Jq — in [32|). 



VI. CONCLUSION 



To avoid the simulation of the dynamics of disordered systems which can be extremely slow, we have proposed 
in this paper to focus on first-passage times that satisfy 'backward master equations'. We have shown that these 
equations satisfy exact renormalization rules upon the iterative elimination of configurations. We have explained the 
similarities and differences with the strong disorder renormalization of Refs [l], ■ We have then tested numerically 
this approach for two types of disordered models : (i) for the random walk in a two-dimensional self-affine random 
potential of Hurst exponent H — 1/2, we have computed the statistics of the first exit time from a square of size 
L x L if one starts at the square center, (ii) for the dynamics of the ferromagnetic Shcrrington-Kirkpatrick model, we 
have studied the statistics of the first passage time tf to zero-magnetization when starting from a fully magnetized 
configuration. We have compared with the results concerning the largest relaxation time towards equilibrium obtained 
with the method of [To| . Our conclusion is that the first-passage method is reliable to measure dynamical properties 
of disordered systems. Although in some cases, it takes more CPU time than the method of [lOj, it can have several 
advantages in other cases : 

(i) it does not require the detailed balance condition (in contrast to 10]) 

(ii) the CPU time depends only on the size of configuration space, but not at all on the disorder realization and 
on the time scales involved that can be arbitrarily large, (in contrast to [l(| where the convergence of the iteration 
method depends on the disorder sample and on the temperature). 

(iii) the freedom in the choice of the initial condition and of the 'target confi gur ations', can be useful to study the 
time scales associated to various dynamical processes (whereas the method of [10] focuses on the largest relaxation 
time towards equilibrium). 
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